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The dimensionally reduced action is believed to provide for a theoretically con- 
sistent and numerically precise effective description of the thermodynamics of the 
quark-gluon plasma, once the temperature is above a few hundred MeV. Although 
dramatically simpler than the original QCD it is, however, still a strongly inter- 
f****) acting, confining theory. In this talk I speculate on whether there could exist a 

£f) i further simplified recipe within that theory, for physically relevant temperatures, 

f**^ ■ which would already lead to a phenomenologically satisfactory description of the 

"""■^ ' free energy and various correlation lengths of hot QCD, but with only a minimal 

amount of numerical non-perturbative input needed. 

Q-t, 



1. Introduction 



It is hoped that one day the reliability of the computational methods de- 
veloped for finite temperature QCD can be systematically checked against 
experimental data from heavy ion collision experiments. Indeed, we could 
then apply the same tools to important problems in cosmology with some 
more confidence. A well-known example of an observable on which a lot 
remains to be understood under both circumstances, is the time scale at 
which a non-Abelian system thermalizes, starting from some complicated 
initial non-equilibrium state. 

We should therefore try to develop methods for determining physical 
observables at high temperatures which are both theoretically consistent, 
numerically precise, and possibly also simple to apply. In this talk, this goal 
is further narrowed in two ways. First of all, we only discuss temperatures T 
sufficiently above the deconfmement transition (or crossover), T>2T ( 9 CD . 
Second, we only consider "static" observables, such as correlation lengths 
and the equation of state. In other respects we do in principle allow for the 
full physical QCD, with Nf > light dynamical quark flavours and a small 
baryon chemical potential [i, although for brevity the data shown below are 
restricted to Nf = 0, fi/T = 0. 



1 



February 1, 2008 20:11 WSPC/Trim Size: 9in x 6in for Proceedings 



proc02'ml 



2 



Since our ultimate interest is not only to be theoretically consistent 
but also numerically precise, we choose for the sake of this talk a rather 
phenomenological approach. That is, for lack of an a priori systematic 
analytic framework for describing confining dynamics, we will compare the 
results of a number of recipes to be introduced presently with actual data, 
although not from collider but from lattice experiments. Thus we obtain 
evidence for some conjectures concerning the non-perturbative dynamics. 

There are also other approaches on the market with somewhat similar 
objectives; for recent reviews see, e.g., Refs. 1 > 2 > 3 . Their general framework 
is however very different from ours, with the exception perhaps of Ref. 4 , 
which also takes the dimensionally reduced theory as a starting point. 

2. Physical scales of the system 
2.1. The asymptotic regime 

Let us start by recalling the parametric scales appearing in finite temper- 
ature QCD in the static limit. The periodicity over the Euclidean time 
direction introduces the scale present already in a non-interacting theory, 
p ~ 2ttT '. For the bosonic Matsubara zero modes which are massless in the 
non-interacting theory, interactions introduce softer scales, corresponding 
to collective plasma modes. In particular, colour-electric fields are screened 
at the Debye scale p ~ gT, while colour-magnetic fields are screened only 
non-perturbatively, at the scale p ~ g 2 T 5,e . Due to three-dimensional con- 
finement, there are no scales softer than this, so that the list is exhaustive. 

As is well-known, perturbation theory has a number of problems in 
this setting. This can perhaps most easily be understood from the fact 
that there is an expansion parameter related to bosonic fluctuations with 
momentum scale p of the form 



where nj is the Bosc-Einstein distribution function. Thus, for hard modes, 
p ~ 2nT, the expansion parameter is ~ 5 2 / 71 " 2 like at T = 0, and should 
remain viable down to a scale 2ttT ~ 2 GeV, which corresponds roughly to 
T ~ 2T ( 9 CD . For the colour-electric modes, p ~ gT, the series is only in 
~ g/n, and seems to converge much worse, unless T 3> T^ cr> 7 ' 8 < 9 > 10 . For 
colour-magnetic modes, p ~ g 2 T, there is no perturbative series 5 ' 6 . 

Thus, there is certainly one simplification from the full theory one can 
make: at T > 2T C QCD , the scales p <~ 2nT can be integrated out. This leads 



-g 2 (2nT)n b (p) = 



g 2 (2irT) p< t g 2 (2irT)T 

■ 



(1) 



7r(eP/ T - 1) np 
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to a dimensionally reduced effective theory for the bosonic zero modes 11 : 

£ E = l - Tr Ft + Tr[A, A a } 2 + m 2 E Tr A 2 + A E (Tr A 2 ,) 2 . (2) 

Here m| ~ g 2 T 2 , Ae ~ g 4 T are generated radiatively 11 . These parameters, 
as well as the effective gauge coupling, have been computed up to next-to- 
leading order (Ref. 10 and references therein). The question we address in 
this talk is whether any further simplifications are possible. 

In the limit of asymptotically high temperatures, corresponding to g <C 
1, the different scales in the system come with the hierarchy 

g 2 T < gT < 2ttT. (3) 

Then the action in Eq. (2) can certainly be simplified, by integrating out Aq. 
One consequence is that all local gauge-invariant operators which do not 
have quantum numbers forbidding this, have correlation lengths of the order 
of the magnetic scale, O^T)' 1 ) 12 . 

2.2. The real world 

For realistic temperatures, however, the situation is different. In fact, the 
longest correlation lengths are associated with operators such as Tr Aq, and 
the corresponding mass scale is significantly below that for purely magnetic 
objects 13 (see also below). Hence the realistic ordering is, in some sense, 

gT < g 2 T < 2ttT. (4) 

It is easy to understand that this has to be the case: close to TjP 00 the 
dynamics starts to be dominated by the Z(iV c ) symmetry related to the 
Polyakov loops, denoted here by Pol (for recent literature, see Ref. 14 and 
references therein), which are represented in the effective theory just by 
Tr A\ — Re Tr Pol, Tr A\ <~ Im Tr Pol. (Note that the discussion here is in 
terms of degrees of freedom represented by spatially local gauge-invariant 
operators, rather than constituents 15 , although this does not really make 
much of a difference) . 

We may thus state, at least on some heuristic level, that: 

• Since the magnetic scale ~ g 2 T is large (cf. Eq. (4) and Fig. 1 
below), it has to be taken into account in all observables where it 
enters additively. 

• At the same time, correlation lengths which are parametrically per- 
turbative (~ gT), remain so non-perturbatively too, because the 
magnetic glueball is too heavy to enter as an intermediate state. 
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3. Conjecture 

Motivated by this physical picture, we propose as an organizing principle 
for the remainder of this presentation the following conjecture: 

If one starts from the dimensionally reduced theory and deter- 
mines physical observables up to and including the leading non- 
perturbative "magnetic" contribution, such that effects from all 
the scales have made their entrances, then one should already have 
accounted for most of the dynamics at the scales ~ gT, g 2 T. 

In the remainder, we discuss numerical evidence in favour of this proposal. 



4. Observables 

The observables to be discussed, and their parametric behaviours, are: 

• The spatial string tension ^fcT s : 

y/o~s ~ [non- pert] x g 2 T+ ... . (5) 

• Correlation lengths £g = related to "eigenstatcs" ~ Tri* 1 ?-: 

TtiQ ~ [non-pert] x g 2 T + ... . (6) 

• Correlation lengths $ = [m^ 1 }^ for ReTrPol, ImTrPol, 
Tr F 0i F jk , etc: 

[m E ] W ~ gT + g 2 T{\n i + [non-pert]) + ... . (7) 

• Pressure, energy density, ... : 

7^7 ~ 1 + g 2 + g 3 + 9 4 In - + .g 5 + ,g 6 (ln - + [non-pert]) + ... . (8) 

t 4 g g 

The conjecture then says that if we determine all the terms shown in 
Eqs. (5)-(8), then the sum of the remaining corrections is maybe on the 
0(20%) level or so, rather than of order unity. 

Let us try to clarify one aspect of the procedure. Indeed, the recipe 
introduced is of course the same as could be used for asymptotically small 
values of g } as in Sec. 2.1. However, in that case the non-perturbative 
correction would be small. For the "real world" case of Sec. 2.2, on the 
other hand, this procedure is just a "recipe" : the expansion is formally not 
sound, since g is not that small and the correction is huge, and it is only our 
conjecture which asserts that the sum of all further terms is subdominant. 
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How are the series shown in Eqs. (5)-(8) to be determined? This de- 
pends on the observable. Formally, one has to systematically integrate 
out A to a sufficient depth, and measure the remaining non-perturbative 
numerical factor, using the theory 

C M = \^Ff j . ( 9 ) 

This is very simple for ^fo~ s and tog, which are obtained from finite observ- 
ables within the theory in Eq. (9), so that one only needs to compute the 
change in the effective coupling constant g. 

For the other observables one in general needs an explicit matching 
computation, because the operators to be used within Eq. (9) are ultraviolet 
divergent. For instance, for electric screening lengths the structure is 12 

[m E ]W = - hm -ln{O^(z)O { i\0))c E (10) 
= : gT + g 2 T In - lim - ln(Q« (z; 0))c m ,a m + 0(g 3 T), (H) 

where Am denotes some ultraviolet regulator within the theory in Eq. (9), 
and Q M ^ (z; 0) is a specific non-local operator, which depends on the corre- 
sponding operator O^(z) in the theory of Eq. (2). The "non-pert" term 
in Eq. (7) arises as a combination of the 2nd and 3rd terms on the right- 
hand-side of Eq. (11). 

Finally, one can carry out a similar matching computation for p(T). 
The first step in the procedure, dimensional reduction, leads to 

£. 1+9Wl »?£ +9 e 1 »?£ + ... + Eg, (12) 



PE 
y4 



= hm JLln [vAiVAoe-f 43 *^. 



(13) 



Here V is the volume, and Ae is the ultraviolet cutoff within the theory 
of Eq. (2). The second part, integration over A 0} leads to 

P E 3 i 4i A E , 5 . 6l A E g TO E p M 

T^= 9 +g ln ^ +9 +9 ln ^ +fl ln A^ + - + Tf (U) 

The terms here were determined up to order 0(g 5 ) in Rcf. 9 , and the 
coefficients of the 0(g 6 ) logarithmic terms are determined in 16 > 17 > 18 . 

Unfortunately, unlike in Eq. (11), the numerical factors inside the log- 
arithms in Eqs. (12), (14) have not been computed yet. This means that 
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the "non-pert" term of Eq. (8) cannot be resolved. Since such numerical 
factors correspond however just to a constant, we can treat this constant 
as a free parameter, and see whether the conjecture has any chance of sur- 
viving. Ultimately, the constant is computable by a collection of various 
techniques, as outlined in Ref. 17 . 

5. Evidence 

We now turn to data. When we discuss the observables in Eqs. (5)-(8), we 
will in general refer to four different kinds of approximations: 

(1) "Full 4d": Results from full 4d lattice Monte Carlo simulations. 

(2) "Full 3d" : Results from 3d lattice Monte Carlo simulations of the 
action in Eq. (2), with the effective parameters expressed in terms 
of the physical 4d parameters and T as specified in Ref. 10 . 

(3) "Leading pcrturbative term" : For expressions of the form 
in Eq. (7), the leading contribution of order gT, given by some 
multiple of the parameter toe of the action in Eq. (2). 

(4) "Up to first non-perturbative term" : Terms up to the order shown 
in Eqs. (5)-(8), in cases where this part has been resolved from the 
"full 3d" result, for instance through a procedure such as the one 
shown in Eq. (11). 

As is well-known, the "full 3d" and "full 4d" results for correlation lengths 
agree almost within statistical errors, and in any case after estimating the 
magnitude of systematic errors arising from lattice artifacts 13 . Our conjec- 
ture then claims that the 4th approximation, which requires only a single 
number of non-perturbative input (determinable with the theory of Eq. (9)) 
rather than full 3d simulations of the more complicated action in Eq. (2), 
should produce an outcome close to the values of the first two procedures. 

5.1. Spatial string tension and correlation lengths 

The numerical data available for the spatial string tension ^fal\ for the 
inverse of the correlation length related to the cigcnstate which connects 
to that determined from TrF?- at asymptotically high T, denoted here 
by m G ; and for inverses of correlation lengths which have a perturbative 
leading term from the scale gT, denoted by [toe]' 1 - 1 ; are shown in Fig. 1. 

We observe that for purely magnetic observables {-^ful, wq), results 
obtained with Eq. (9) agree nicely with full 3d results, as well as with 
full 4d results, where available. In particular, even though in principle 
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10 x (a )"V 

m G /T,TrF |2 F |2 , SU(2) - 

m G /T,TrF |2 F |2 ,SU(3) - 

m/r,TrF 0| F 23 ,SU(2) - 

m E /T,TrF (|| F 23 ,SU(3) - 

m E /T,ReTrPol,SU(2) - 

m/T, ReTrPol, SU(3) - 

m/T, ImTrPol, SU(3) - 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 

Figure 1. A comparison of spatial string tensions, as well as inverses of correla- 
tion lengths in various channels and two gauge groups, at T m 2T C QCD , with dif- 
ferent methods, as explained in the text. The lattice data have been collected from 
Refs. !3, 19, 20, 21, 22 Representative continuum operators have also been shown. 

there, there is no need to include any higher order corrections beyond those 
determined by Eq. (9). This observation is certainly in perfect agreement 
with our conjecture. 

For the electric observables, on the other hand, the "up to first non- 
pert, term" -prediction has so far been resolved in only one quantum number 
channel, although it could in principle be resolved in others, as well. We no 
longer observe perfect agreement. However, while the leading perturbative 
terms are too small by a factor ~ 3, "up to first non-pert, term" predictions 
arc too large by ~ 30 %. Thus the conjecture still works semi-quantita- 
tively: all the higher order corrections in Eq. (11) sum up to a small number, 
even though the first correction beyond the leading term is huge. It would 
be interesting to see whether the conclusion holds also after studying other 
channels, where the non-perturbative correction should be much smaller, 
in order to fit the full 4d and full 3d results (cf. Fig. 1). 

5.2. Pressure 

We finally turn to the pressure. Unfortunately the situation remains some- 
what open here. Not only has the non-perturbative g 6 term in Eq. (8) not 
been resolved, but we also do not have a "full 3d" prediction available: 
measurements can be carried out but a constant again remains open 24 > 25 . 
These problems are ultimately due to the fact that the constant parts re- 
lated to the divergent 4-loop matching coefficients in Eqs. (12), (14) have 



H2H 
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□ full 3d 

[> leading pert, term 

< up to 1st non-pert, term 
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1 10 100 1000 

T/A MS 

Figure 2. The pressure, normalised to the non-interacting Stefan-Boltzmann law po, 
for various values of the 0(g 6 ) constant term (from Ref. 17 ; see there also for further 
details on what exactly is plotted here). The 4d lattice data is from Refs. 19 > 23 . 

not been determined yet. The only unambiguous results are from 4d lat- 
tice simulations (at least for Nf = 0), and from perturbation theory at 
asymptotically high temperatures. 

Therefore, the best we can do is to see whether a necessary condition 
for the validity of the conjecture is satisfied: is it possible to find some 
value representing the unknown 4-loop constant terms, such that it repro- 
duces a whole function of the expected type? This exercise is carried out 
in Fig. 2 (taken from Ref. 17 ), and we indeed find that one can reproduce, 
in principle, a reasonable curve with a value of the constant of order unity. 

Of course, it must be said that due to the high perturbative order ap- 
pearing, the pressure is an observable which is quite sensitive not only to 
the non-perturbative term but also to ambiguities in perturbation theory. 
For illustration, we demonstrate this in Fig. 3. An optimal value has been 
chosen for the constant, and then the renormalisation scale entering the 
expressions of the effective parameters in Eq. (2) has been varied within 
the range 0.5/2 O pt...2.0^ O pt, where fL opt ~ 6.7T is suggested by the next-to- 
leading order expression for c/| 10 . 

The scale dependence is non-monotonic, because the first derivative of 
the result with respect to p, is very close to vanishing, due to the fact 
that we have gone via the effective theory in Eq. (2) to obtain the output. 
Nevertheless, we see that at T < 3T C QCD (T C QCD ~ A m ) the predictions are 
in any case getting unreliable, irrespective of the value of the constant. This 
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Figure 3. The dependence of the result of Fig. 2 on the renormalisation scale /i, around 
Aopt ~ 6.7T, for an 0(g 6 ) constant ~ 0.7 (cf. Rcf. 17 ). 

may not be too surprising, though, since to determine the full g 6 constant 
term requires also renormalisation of g at 2-loop level, but this has not been 
incorporated in the expressions shown here. 

6. Conclusions 

There is some numerical evidence for the conjecture of Sec. 3 from the 
spatial string tension as well as from various ("magnetic" and "electric") 
correlation lengths, but more channels should be added. It would also be 
interesting to consider other observables (see, e.g., Ref. 26 ). 

The conjecture might also work for the pressure, although so far only 
one necessary condition is observed to be satisfied, and a lot of work remains 
to be carried out before we know whether a sufficient condition is there. 

Once such analytic expressions with non-perturbative constants are 
known, extensions to finite Nf > and n<irT are easy (see, e.g., Rcf. 27 ). 

Hopefully similar expressions can be worked out also for real time ob- 
servables, where 4d lattice simulations are even much more difficult. 
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